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We extend the composite fermion construction to the torus geometry. We verify the validity of our 
construction by computing the overlap of the composite fermion state to the exact diagonalization 
ground state of both Coulomb interaction and Haldane-pseudopotential interaction Vo (Vi) for 
bosonic (fermionic) states. 
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I. INTRODUCTION 

A paramount impetus for the growing interest in 
strongly correlated quantum matter is the discovery that 
such systems can be topologically ordered. The first, and 
most prominent, examples are the various incompressible 
fractional quantum Hall (FQH) liquids. These are formed 
at low temperatures, when very clean two-dimensional 
electron gases are subjected to a strong perpendicular 
magnetic field [1J. The most striking consequence of the 
topological order in the FQH liquids is that the emergent 
low-energy quasiparticles have fractional electric charge, 
and are believed to obey fractional braiding statistics 
[5] • Another hallmark of a topological phase of matter is 
that the number of degenerate ground states depend on 
the topology of the space on which the state is denned. 
Again, the simplest example are the FQH liquids, where 
numerical solutions of the microscopic Coulomb problem 
can be compared to different theoretical predictions. In 
this paper we develop techniques that makes it possible 
to make such comparisons for the Jain states, which are 
an important set of actually observed FQH liquids. 

The problem of many particles moving in a strong 
magnetic field, and interacting via Coulomb forces, is 
intractable, so one has to find an effective description. 
While a natural approach is to construct effective low- 
energy field theories |3_ , a very fruitful alternative route 
— following a seminal article by Laughlin [3] — has been 
to construct model wave functions with well defined 
topological properties, and verify numerically that they 
describe small systems accurately. The most common 
method of verification is to compute the overlap with the 
exact Coulomb ground state, but recently it was shown 
that studying the entanglement spectrum |S] can give 
valuable complementary insights. 

Although most effort in FQH physics during the last 
decade has been aimed at understanding different non- 
abelian states, we shall here concentrate on the family 
of observed states in the lowest Landau level. There are 
two successful theoretical approaches to describe these 
states — the Haldane-Halperin hierarchy (6j |7] and the 
composite fermion (CF) theory [5]. 

The Haldane-Halperin hierarchy describes a family of 
incompressible FQH liquids formed via successive con- 
densation of the low-energy excitations — quasiholes and 



quasielectrons. It predicts that incompressible FQH liq- 
uids may be found at filling fractions v = p/q, where p 
and q are relatively prime integers and q is odd. It also ar- 
gues that the stability of the liquids decreases roughly as 
~ 1/q with increasing denominator. The emergent quasi- 
particles have fractional electric charge 1/q and obey 
fractional abelian braiding statistics. The simplest way 
to obtain explicit wave functions for ground and quasi- 
particle states, at all levels of the hierarchy, is by using 
conformal field theory techniques [HI HO] • 

The CF theory can describe most of the observed FQH 
liquids by mapping the problem of strongly interacting 
fcrmions that fill a fraction v of a Landau level, to that of 
non-interacting (or at least weakly-interacting) compos- 
ite fermions filling an integer number of effective Landau 
levels in a reduced magnetic field. The latter describes an 
incompressible state because of the finite gap between the 
effective Landau levels. This approach also gives a simple 
picture of the ground state and the low-energy sector. A 
big initial success of the CF theory was the very good 
agreement with the exact Coulomb eigenstates obtained 
from finite-size numerical studies, both for the ground 
state and the excited states. It was later shown that for 
the important case of the positive Jain series, the hierar- 
chy wave functions exactly coincide with those from CF 
theory, both in planar and spherical geometry |11H13| . 
The two approaches also give very similar predictions for 
the relative stability of the various FQH liquids. 

Both approaches to the abelian FQH states in the low- 
est Landau level — the hierarchy and the CF theory — 
are well-understood (and studied) in the disk and sphere 
geometry, see eg. |14j and references therein. The torus 
geometry, however, has been studied far less and explicit 
wave function were only known for certain model states, 
that are determined uniquely (up to center-of-mass trans- 
lations) by their vanishing properties [TBl - fTT] . A first 
attempt to construct hierarchical wave functions on the 
torus was done in Ref. |18j using conformal field theory 
techniques. While it provided wave functions that are 
very good approximations to the exact Coulomb ground 
state in a certain parameter regime, the construction was 
not satisfactory in that the wave functions did not trans- 
form properly under modular transformations and were 
uniquely defined only in the thermodynamic limit. It 
has only recently been understood how to resolve these 
problems |19j . 



In this paper, we show how to generalize the CF con- 
struction to the torus geometry. The construction has no 
free parameters and gives unique model wave functions 
(up to center-of-mass translations) of the ground state 



and the excited states at filling fraction v = 



np+l ' 



and 



p integers, corresponding to the positive Jain series. 

There are several reasons, why the torus geometry is 
interesting, even though it cannot be realized experi- 
mentally. We already mentioned the topological ground 
state degeneracy, but it is also important that numeri- 
cal calculations are better defined on closed manifolds, 
such as the sphere and the torus, since they do not suf- 
fer from edge effects, which can be substantial for the 
system sizes one can reach numerically. While numer- 
ics on the sphere have proven very useful, there are still 
problems connected to finite size, most notably the so- 
called shift. It can happen that two states that are at the 
same (thermodynamic) filling fraction v appear at differ- 
ent magnetic fluxes jV^ in the finite-size system. The 
most prominent examples are the Moore-Read state [50] 
and the CF Fermi liquid [5T], which are both at filling 
1/2 but have different shift. Numerical comparison of 
these two states on the sphere is therefore only indirect. 
On the torus, this issue does not arise, which allows for 
a directly comparison of these states [22 . 

Another advantage of the torus geometry is that one 
can change the shape of the torus — described by the 
modular parameter r — and thus get more information 
about a state without having to increase the system size. 
This was successfully used for entanglement entropy cal- 
culations, where one wants to extract a subleading con- 
stant term in the entropy. Ref. [53J showed that changing 
the aspect ratio of the torus and thus obtaining additional 
data, yielded much more accurate bounds on the topolog- 
ical constant than could be obtained from sphere calcu- 
lations. Also, as shown by Avron et.al. [53], by studying 
the response of QH liquid to an adiabatic change in r, one 
can determine the odd part of the viscosity tensor. An 
explicit calculation in the case of the Laughlin states was 
made by Read [55] , and the results in this paper could be 
used to perform similar calculations for the Jain states. 
Let us also note that the techniques introduced in this 
manuscript are not restriced to the positive Jain series, 
but can also be used to study more exotic states, such 
as the BS states [55], the NAC states [57], as well as the 
closely related bipartite CF states [55] . 

Outline of the manuscript In Section ITT] we first 
present the CF construction on the disk geometry and 
show how to generalize the approach to the torus geom- 
etry. We discuss single-particle states on the torus in 
Subsection |II B| and give an expression for the product of 
two such states at different magnetic fluxes N^ and N^, 2 
in Subsectio n |IIC| The derivation of this identity is given 
in Appendix |A| In Subsection |II D| we show how to eval- 
uate CF states on the torus. Overlaps of some CF states 
with exact diagonalization results using both Coulomb 
and Haldane pseudo-potential interactions are calculated 
in Section |III| These overlaps should be regarded solely 



as a proof of principle that the construction on the torus 
is sound. In Section |IV| we speculate on possible lowest 
Landau level projection schemes in real space, given that 
the torus places additional constraints on model wave 
functions. 



II. GENERAL COMPOSITE FERMION 
CONSTRUCTION 

In this Section we explain how to generalize the CF 
construction to the torus geometry. In Subsection II A we 



first discuss the CF construction on the disk and sphere 
and point out some subtleties that become important 
on the torus. Subsection IIIBI contains a short review 
on the single-particle states on the torus. In Subsection 
|HC| we derive formulas for the projection of a product of 
two single-particle states. In Subsection |IID[ we discuss 
properties of the CF states on the torus using the bosonic 
CF state at v = 2/3 as an explicit example. 



A. Composite fermions on the disk geometry 

There are already many good texts on the CF 
construction — see for instance [T3] for an extensive and 
pedagogical review. Thus, we keep the discussion in this 
section very brief and focus on properties that are rel- 
evant for the torus. In the following, we restrict our- 
selves to the positive Jain series at fillings v = " 1 , 
where n > 1 and p are integers. We expect the negative 
Jain series to work analogously, but we have not yet per- 
formed any explicit calculations. In the CF theory, one 
attaches an even (odd) number of vortices to strongly in- 
teracting fermions (bosons). The resulting fermionic par- 
ticles are called composite fermions and one assumes that 
these composite particles are non-interacting or at least 
very weakly interacting. Due to the attachment of vor- 
tices, they feel a reduced magnetic field B* = B — ppcfio, 
where </>o is the magnetic flux quantum and p is the two- 
dimensional density. For properly chosen p the reduction 
of magnetic flux is such that the CFs fill an integer num- 
ber of effective Landau levels. 

A trial wave function for the ground state of strongly 
interacting particles at filling v — " 1 is then usually 
written as 



M{*i» = VLLL I *n({*i» • IIC* - ztf \ , (1) 

where z — x + iy is a complex coordinate. &n({zj}) 
is the many-body wave function (slater determinant) for 
the n lowest Landau levels filled and Vlll projects to the 
lowest Landau level. Eq. §\l does not strictly speaking 
describe a proper lowest Landau level wave function on 
the disk, because it does not have the correct gaussian 
factor. Usually, one does not worry about this but just 
adds the correct factor by hand. However, this subtlety 



becomes important on the torus as explained in the next 
paragraph. 

The naive guess of how to generalize Eq. (fTj), is to 
replace each part by the respective torus counterpart. In 
particular, this would amount to replacing the Jastrow 
factor with its periodized version [15] 
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where d\ is the odd Jacobi theta function (defined by set- 
ting a = b = 1/2 in Eq. ( 15 1). We choose Q\ because it is 
the only ^-function that has the correct short-distance be- 
haviour, ie. is it the only antisymmetric 6- function. How- 
ever, this choice poses two obvious problems: firstly, the 
wave function does not obey the correct boundary condi- 
tions on the torus, see Eq. dTl. [37J Secondly, there is no 
efficient way to project the many-body wave function to 
the lowest Landau level. To the best of our knowledge, 
no analog of the Girvin-Jach projection [29] is known on 
the torus. We will comment more on this in section HVl 
Instead of Eq. (fTl) we will consider the following ex- 
pression: 



*„({*i» = V LLL {**({*,}) • <&!({*,-})*} 



(3) 



The replacement of the Jastrow factor to $i is of course 
trivial for both the disk and sphere — in the former case 
it only differs by a gaussian factor. The point is, how- 
ever, that ([3j) is a, proper Landau level wave function on 
the disk, ie. it has the correct gaussian factor because 
the gaussian factors of <& n and Q\ combine to give the 
correct factor at the combined flux. In addition, using ex- 
pression fcfy solves both the problems mentioned in the 
previous paragraph. It is straightforward to verify that 
&u({zj}) obeys the boundary conditions on the torus d7|. 
The projection onto the lowest Landau level can be im- 
plemented on the single-particle level, which is explained 
in subsection III CI 



B. Single-particle states on the torus 

We consider a torus spanned by two, not necessarily 
orthogonal, translation vectors L\ and L%. A homo- 
geneous external magnetic field — perpendicular to the 
surface of the torus — is described in terms of the vector 
potential A = —Byx using Landau gauge. The number 
of flux quanta piercing the torus is related to the area 
A= |Li xL 2 \ of the torus by 2-Kt B N (j) = A = L X L 2 sin(6>) 
with magnetic length £ B = \Jhc/{eB) and 9 being the 
angle between L\ and L 2 . The case of a rectangular 
torus corresponds to — it/2. The shape of the torus is 
conveniently parametrized by the aspect ratio 



v v* . 



Lx 



(4) 



In the presence of the magnetic field, any valid wave 
function on the torus must be invariant (up to an overall 



phase) under single-particle magnetic translations t(L{) 
and t(L 2 ), where the magnetic translation operator is 
defined as 



i(L)=explL-(V~i^A)-t L ^-}. (5) 



Let us define 'small' magnetic translations 




x exp 



(6) 



The periodic boundary conditions of a wave function ip 
can, thus, be formulated as: 



i N <t> i 
t N -t> i 



2 1p. 



(7) 



In the remainder of the manuscript, we will set the 
solenoid fluxes a\, a 2 — without loss of generality. 

As t\ and t 2 do not commute with each other, we can 
choose the single-particle states to be eigenstates of only 
one of them. In the following, we will mostly use eigen- 
functions of t\ 



ifa{x, y )=K B E e_2 " 



fc= — oo 



x exp 



N r , 



-U + kNtf 






(8) 



where z = (x + iy)/L± is the dimensionless complex co- 
ordinate of the particles, n = 0, 1, ... is the Landau level 
index and j = 0, . . . , (A^ — 1) the momentum index. Note 
that the momentum is only defined modulo N^, because 
any larger value can be absorbed into the sum over wind- 
ings around the torus. The normalization constant is 
given by 






y/2N^(r) 
(2™nW) 



1/2 



where A 



2irt B Nt 



is the total area of the torus. 



(9) 



n,j 



is an eigenfunction of t\ with eigenvalue exp[— 2nij /N<p], 
while t 2 shifts the momentum by one: t 2 <f> r f- = <\>^j-\- 
As we will need to distinguish single-particle states at 
different flux later on, we keep the magnetic length i B as 
an explicit parameter in the single-particle state </>„ s 7 . 



C. Product of single-particle states on the torus 

In complete analogy to the disk and sphere geometry, 
we can write a product of two single-particle states on 
the torus at magnetic flux N ( j )1 and N<f, 2 as 



I11+T12 N^ — l 

E E 

n=0 j=a 



CM"-4>.¥). (10) 



where the magnetic lengths are related by 



(: 



2 ! 



which is equivalent to N,h 



N, 



N* 



The con- 



stants Cl)%j n = Cl)^l n (N^N^,r) depend on the 
fluxes Nfa 2 as well as the aspect ratio of the torus. They 
can be computed for arbitrary n\ and n 2 , but we have 
not been able to find a closed formula except for very 
simple cases. For the CF construction we need to know 
the coefficients for n 2 = n = 0, but arbitrary ri\. In the 
following, we restrict ourself to these cases. 



In order to simplify notation later on we define Q as 
the greatest common divisor (gcd) of JV^ and Nq 2 : 

Q=gcd{N 4>1 ,N 4>2 ) 
N^ =h-Q 
Nfc = t 2 -Q 
N <l> = (t 1 +t 2 )-Q = t-Q. (11) 

It follows that Q = gcd(N^ >1 ,N (t> ) = gcd{N <j>2 ,N (j> ). The 
different magnetic lengths are related to t\ and t 2 by: 



(12) 



For ni =0,1 the coefficients in ( 10 1 become rather 
simple: 
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(13) 
(14) 



where j = {ji+j 2 +(3Qti) mod JV^, for fl = 0, 1, . . . ,t— 1. 

<%&*r c hSu = ° for ** 0*1 + h + fiQti) mod n*. 

The third ^-function is defined by: 



»(*!«) 



z\q) 

V^ Jk+a) 2 e 2i(k+a)(z+b) 



(15) 



and 9' 3 (z\q) = d z 9 3 (z\q). For higher n 1; the coefficients 
can in principle still be represented with help of higher 
derivatives of the ^-function, but they become increas- 
ingly cumbersome to evaluate. They can be written as: 
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(16) 



The derivation of Eq. (16 1 involves straightforward, but and (14) can be obtained from (16 1 by a Poisson resum- 
tedious algebra, which is done in Appendix [A! Eq.s (13) mation. Note that for 3(r) not too small, the summa- 



tion over s converges rapidly. For numerical purposes, 
one needs to care only about the first few terms around 
zero. The problem of convergence for small 9(t) can in 
principle be avoided by doing a Poisson resummation on 
the sum over s — similar to what was done in obtaining 



Eq.s (131 and (14 1 



D. Composite fermion wave functions on the torus 



describes a fermionic (bosonic) state. As a sanity check, 
one may note that n = 1 reproduces the Laughlin states 



at filling -^pj. Expression (|T7|) can also be used to eval- 



uate wave functions corresponding to quasihole and/or 
quasielectron excitations — in the same way as on the 
sphere or the disk. For sake of simplicity, we focus on 
ground state wave functions in the following discussion. 



The formulas derived in the previous section allow for 
evaluation of the Jain state at filling v = — ^r given by 



*.({^})='PiLL [*„({*,-}) •*l({« i }) P ] 



(17) 



where $., is the many-body wave function of the lowest 
j Landau levels completely filled. For p even(odd), this 



In principle, it is straightforward to evaluate (17 1, by 



multiplying out the slater determinants and using the 
coefficients ( [T6| (repeatedly if p > 1) to reduce the ex- 
pression to a lowest Landau level wave function at the 
combined flux. For the simplest state, describing the 
bosonic Jain state at filling v = 2/3, the explicit expres- 
sion becomes: 



' N-l 2 



* a/ 8({*i,*})= e(-d ct n e n (cu^c*. 



<tGS a 



+N/2);j a+N/2 
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M J] njlm^izj}), (18) 



where j a — (a modA^ + cr(a) + f$ a Qt\) mod N^. The 
rn fj.({ z j})' s are the many-body basis states on the torus: 



m M z j}) = (zi, . . . , z N \n) , 



(19) 



with /j, = {n , . . . , riN^-i} and n 3 = J2a=o $jj a being 
the occupation number of the single-particle orbital with 
^-eigenvalue exp[— 2nij/N^\. 



The qualitative difference between Eq. ( 18 1 to the cor- 



responding disk and sphere expressions lies in the addi- 
tional sums over j3i , . . . , /3jv- In the disk and sphere geom- 
etry, the momentum of the product of two single-particle 
states on the torus is simply the sum of the two momenta. 
The projection, thus, involves evaluating N\ terms in or- 
der to obtain the coefficients in the occupation number 
basis. On the torus, the momentum is only defined mod- 
ulo the flux Nfj,. This implies that the winding sums in 
the single-particle states at flux N ( f )1 and N<p 2 yield dif- 
ferent momenta at the final flux. For instance, in order 
to compute the CF state at filling fraction v = 2/3, one 
needs to evaluate ^13^ terms, which limits the system 
sizes one can reach. Note that this limitation becomes 
worse, if we increase p in Eq. (17). 



In addition, ( 17 1 is in general not an eigenstate of the 



many-body translation operator T\ = Ylj=a n > wnere 

if translates the j'th particle by Lx/N<f,. An eigenstate 
can be obtained by either restricting to the correct mo- 
mentum sector in the Fock basis or by applying the ap- 
propriate projection operator. 




Figure 1: color online: Overlap of CF state at filling v — 2/3 
for N — 10 bosons with the exact diagonalization state using 
Coulomb interaction (red x ) and a contact interaction using 
the Haldane pseudopotential Vo (blue +). The lines are a 
guide to the eye. 



III. NUMERICAL ANALYSIS 



In this section, we present numerical checks on the 
wave functions obtained by Eq. (17). We computed the 



overlap between the exact diagonalization ground state 
and the bosonic CF state at filling v — 2/3 for system 
sizes up to 10 particles and in the fermionic case at filling 
v = 2/5 for system sizes up to N — 6 particles. The exact 
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tion of these dips depends on the system size — in the 
bosonic case, it seems to be shifted to lower values of \t\ 
for increasing system size. In the fermionic case, there 
is too little numerical data to make a statement. The 
origin of the dips is not clear at the moment, but one can 
note that the overlap of the fermionic Laughlin state at 
v = 1/3 with the Coulomb ground state has a qualita- 
tively similar behaviour as a function of \t\ to the one 
shown in figure [2] 



IV. PROJECTION SCHEMES IN REAL SPACE 



Figure 2: color online: Overlap of CF state at filling v = 2/5 
for N = 6 electrons with the exact diagonalization ground 
state using Coulomb interaction (red x) and a short-range 
interaction using the Haldane pseudopotential Vi (blue +). 
The lines are a guide to the eye. 



diagonalization was done both for Coulomb interaction 
and the smallest relevant Haldane pseudopotential |6J— 
Vo for bosons and V\ for fermions. The shape of the torus 
was kept rectangular (8 = -| in Eq. Q), with aspect ra- 
tios \t\ varying from 0.1 to 10. Due to the invariance of 
the shape of the torus under the modular transformation 
S : t — > — 1/t, we can restrict the analysis to aspect ra- 
tios 0.1 < \t\ < 1 without loss of generality, as a torus 
with \t\ > 1 can be obtained by an ^-transformation. 
The absolute value of the overlap O = K^cFl^ex)! be- 
tween the CF state and the exact diagonalization ground 
state is shown in figures [l] and [2] for varying aspect ratios. 

When comparing overlaps of the torus and sphere ge- 
ometry, we choose the most isotropic point, namely the 
square torus with |r| = 1. In the bosonic case, we 
find an overlap of O — 0.996 (Coulomb interaction) and 
O = 0.973 (Vb-interaction) for a square torus, which are 
slightly higher than the overlaps found in Ref. [30 for the 
spherical geometry. In the fermionic case, we find over- 
laps O = 0.999 (Coulomb interaction) and O = 0.990 
( Vi-interaction) , which are slightly lower, but still com- 
parable to the overlaps found in Ref. [3T] . 

The overlaps depend strongly on the shape of the torus, 
even though they remain quite high throughout the whole 
range of aspect ratios. On general grounds, we expect 
the overlap to approach unity in the limit of \t\ — > and 
\t\ — >• oo. In Ref. [22] it was shown that the ground state 
of Hamiltonians with quite generic repulsive interactions 
becomes a product state for aspect ratios \t\ — > oo — the 
so-called thin torus limit. We checked numerically that 
the CF state has this property as well. The other limit 
\t\ — > can — for a rectangular torus — be mapped to the 
thin torus limit, when using the Landau gauge A = Bxy 
and eigenfunctions of the ^-operator. 

For all system sizes, we observer a dip (in case of 
fermions several dips) in the overlap curve. The posi- 



In principle, the method described in this manuscript 
allows one to compute any CF states on the torus. How- 



ever, evaluating Eq. ( 17 1 becomes numerically hard for 
large p when using Eq. (10 1. On the torus, one needs 
to evaluate (Nl) p t N terms to obtain the wave function 
in Fock space. Even on the disk and the sphere, where 
one does not have the additional complication of the 
/?!,..., (3n sums, the number of terms, which one needs 
to evaluate, is still (Nl) p . This restricts the system size 
to very small systems for large p. 

A way around this, at least on the disk and the sphere, 
is to evaluate Eq. ( 17 1 in real space and use Monte Carlo 



techniques to study the resulting model wave functions. 
How to write the projection operator Vlll in real space 
was shown by Girvin and Jach in Ref. [29 . It amounts to 
moving all anti-holomorphic components to the left and 
replacing them by derivative operators z — > 24- with 
the assumption that the derivatives do not act on the 
Gaussian factors. 

This simple implementation of the projection has no 
straightforward analog on the torus. Note that deriva- 
tives are not valid operators on the torus, because they 
destroy the periodic boundary conditions (l7|. Follow- 
ing Ref. [18J one may argue that the torus analog to 
the derivative operator should be related to the small 
translation operators i\ and £2 (Rj), because the small 
translation operators keep the boundary conditions in- 
tact and become effectively holomorphic derivatives in 
the limit Li,L2,N c / 1 — > 00. In fact, the analog of the 
derivative operators may be given by a sum of transla- 
tion operators that have good modular properties [19J. 
It is an open question whether or not this idea could 
be used also to simplify the evaluation of the CF torus 
wave functions. As a modular invariant sum necessarily 
involves iV? terms [33 , this projection scheme would be 
numerically cheaper than the one used here only if the 
sum over different translations converge rapidly. 

A numerically very efficient way to evaluate Eq. 
approximately is the Jain-Kamilla projection 
The Jain-Kamilla projection is a close approximation to 
the exact projection, but can — in contrast to the latter — 
be evaluated for very large system sizes. It amounts to 



dividing the Jastrow factor in Eq. (TTT) as 
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and multiplying each column a of the slater determinant 
*&n({zj}) by EL^a^a — z j) P ■ Then, each component 
of the slater determinant is projected separately to the 
lowest Landau level by 
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(21) 



When trying to generalize this scheme to the torus, one 
may note that the effect of the derivative operators in 



(21 1 is to remove zeroes between particles. On the torus, 



we cannot remove zeroes but we may shift them. This is 
in fact exactly, what the translation operators do, that 
were mentioned in the previous section,. However, one 
needs to shift the zeroes without destroying the boundary 
conditions. The most straightforward implementation of 
(21 1 would, thus, require a doubly periodic function with 



only a single pole, which — as we know from complex 
analysis — does not exist. Unfortunately, the boundary 
conditions impose rather strict rules on how one may 
change the wave function, which is why we have not been 
able to find an analog of the Jain-Kamilla projection on 
the torus. 



SUMMARY 



comparable to the ones obtained in the disk and sphere 
geometry. It turns out that numerical evaluation of the 
wave function is harder on the torus than on the disk and 
sphere, because the winding sums mix different momen- 
tum sectors. We have also speculated on possible gen- 
eralizations of the real space projection schemes to the 
torus geometry, though, unfortunately, we have not been 
able to find an explicit realization. Such schemes may 
allow one to reach larger system sizes than are possible 
with the method presented here. 

Let us emphasize again that our method works for the 
whole low-energy sector of the CF states, even though we 
only treated the ground states explicitly in this article. 
The techniques introduced here may be useful for systems 
that cannot be studied directly on the sphere, because 
they have different shifts, as eg. the one studied in Ref. 
|36] . They can also be used to study generalizations of 
the abelian Haldane-Halperin hierarchy, such as the BS 
states [26 , the NAC state [27 , or the bipartite CF states 
[28] . Also it will clearly be interesting to see whether the 
exact agreement between the hierarchy and the CF wave 
functions that have been demonstrated on the plane and 
on the sphere, also holds true on the torus. 
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In this manuscript, we generalized the CF theory to the 
torus geometry. We showed the validity of our method 
by calculating the overlap between the CF states and 
the exact diagonalization ground state of Coulomb and 
the smallest relevant Haldane-pseudopotential interac- 
tions for filling fractions v = 2/3 and v = 2/5 and 
system sizes up to N — 10 (6) particles for the bosonic 
(fermionic) states. The overlaps on the square torus are 



Appendix A: Derivation of product formula 



In this appendix, we derive formula (|16| for an expan- 
sion of the product of two single-particle states at fluxes 
N ( p 1 and N^, 2 in terms of the single-particle states at the 
combined flux N^ with N^ >1 + Nfc = N^. The product 
of two single-particle states is given by: 



CO oo 

k±— — oo k 2 — — oo 

^2 i a_ , u. at ^2 



x exp 



VKT 



Ui + hN^y {j2 + k 2 N 4l2 y 



N rh 



N, 



<f>2 



■H ni {^( Jl + k 1 N, 1 )- ] L 



(Al) 



In the following, subscript 1 and 2 denote quantities of cate things unnecessarily. 

the two single-particle states respectively, while those 

without subscript denote those of the product. Note that 

Nfj, — N,p 1 + Nfo and £~ 2 — f^ 1 + £% . Let us assume 

n2 = but m may be arbitrary for the time being. We 

could also consider arbitrary n 2 , but it will only compli- 



Define Q as the greatest common divisor (gcd) of N c / >1 
and JV«i„: 



Let us now consider the factor that is exponential in 



the winding number. Using (A7l it can be rewritten as 



Q=gcd(N^,N^ 2 ) 

Nfa = t 2 -Q 

N < f, = (ti + t 2 )-Q 



t-Q. 



(A2) 



It follows that Q = gcd(N 4>1 ,N 4> ) = gcd{N$ 2 , N4). We 
use that the different magnetic lengths are related as : 



h 
t 
T 2 



(A3) 



We define the torus as in Section |TTJ It is easy to check 
that the product of the two single-particle states obeys 
the correct boundary conditions for flux N$. 

Let us first discuss how to rewrite the double sum 
over windings coming from both single-particle states, 
denoted by k\ and k 2 . We can choose integers k, s 6 N, 
and j3 <E {0, . . . , t — 1} such that 



which implies 



k\ — f3 + k — t 2 S 
k 2 = k + hs, 



t-k = t 1 k 1 +t 2 k 2 + /3 
t ■ s — j3 + k 2 — k\ . 



(A4) 



(A5) 



It is beneficial to introduce some more notation, that will 
simplify expressions later on. We find that we can rewrite 
the phase factors of the single-particle states as 



ji + fciiVfc = -^A k - ht 2 Q ■ Y s 
n 

£ 2 
j 2 + k 2 N^ = -^A k +t 1 t 2 Q-Y s , (A6) 

where A k depends only on k and /3, but not on s, while 
Y s depends only on s and /?, but not on k: 

A k =]i+J2 + kN 4> +pt 1 Q. 
hji - hj 2 + Pt\t 2 Q 



Y* 



htzNt 



(A7) 



Using these definitions we see that the z-dependent factor 
on the right-hand-side of Eq. (Al) does not depend on 



the summation index s and becomes rather simple: 

e -2TTi(j 1 +k 1 N 4>1 +j 2 + k 2 N^ 2 )z _ e ~2iriA k z (A8) 



exp 



Ui + kiN^) 2 {j 2 + k 2 N^f 



N d 



exp 



N^ 2 



(A9) 



ie. it factorizes into two parts, each of which only depends 
on one of the summation indices. 

Most of the complication lies in the Hermite polyno- 
mials, at least if n\ 7^ 0, 1. The Hermite polynomial in 
Eq. ( Al ) can be written as 
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\ L x 
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With the following identities 



H n (x + y) = J2{ k )Hk(x)-(2y) 

Lfe/2J 

Hkdx) = £ 7 fc - 2 *( 7 2 - 1) 

i=0 



we find that 
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xi ■iuui'"-'.//, -'(xr^l !AI2) 



where a\ = 2ir£it 1 t 2 Q / L\ and \kj2\ — k/2 if k is even, 
resp. (k - l)/2 if k is odd. The first line in Eq. (TaTTi 



can easily be derived using the generating function of 
the Hermite polynomials. The second is more tedious to 
derive and follows from successive partial integration of 
the integral J dxexp[— x 2 ]H n ( r )x)H m (x). 



niO;n 



We can now identify the coefficient C- 1 ■' . ■ to be 
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(A13) 



where j = (j 1 + j 2 + /3tiQ)m.od N ^. The coefficients setting n = 0. In order to find Eq.s (13 1 and (14 1 one 
vanish for other values of j. Eq. (16 1 is obtained by must do a Poisson resummation on the sum over s. 
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